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Exceptionally bright quasars with redshifts up to z=6.28 have recently 
been discovered^. Quasars are thought to be powered by the accretion of gas 
onto supermassive black holes at the centres of galaxies. Their maximum 
(Eddington) luminosity is proportional to the mass of the black hole, and so 
these bright quasars are inferred to have black holes with masses of more than 
a few billion solar masses. The existence of such massive black holes poses a 
challenge to models for the formation of structures in the early Universe, as 
it requires that the black holes would grow so massive in less than a billion 
years after the Big Bang. Here we show that up to a third of known quasars 
with z~6 will have their observed flux magnifled by a factor of 10 or more 
through gravitational lensing by galaxies along the line of sight. The inferred 
abundance of quasar host galaxies, as well as the luminosity density provided 
by the quasars, are therefore substantially overestimated. 

The four highest redshift quasars known^ with z ^ 5.8 (SDSS 1044-0125 was later 
found^ to have z = 5.73), were selected in the SDSS photometric system to have magnitudes 
z* < 20.2 and colors i* — z* > 2.2. The masses of the black holes powering these quasars 
are estimated to be 3 x IO^Mq, implying^ that the mass of their host galaxies is 
^ 10^^ Mq. Such massive hosts lie on the steep exponential tail of the Press-Schechter^ 
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mass function, so that a correction to the inferred black-hole mass severely impacts the 
estimated cosmological density of galaxy halos which are sufficiently massive to host the 
observed quasars. 

Gravitational lensing leads to a magnification by a factor /j, of the apparent source 
luminosity. Let us first consider the implications of this magnification for the inferred 
properties of the massive quasar systems within structure formation theory. Using the 
Eddington luminosity to set a lower limit on the inferred black hole mass^, one finds that 
the inclusion of the magnification due to lensing lowers the minimum black hole mass by 
a factor of /i. This, in turn, implies that the black hole can form in a lower mass galaxy. 
Figure 1 shows the resulting enhancement factor in the space density of galaxy halos that 
could host the observed SDSS quasars at 2; = 6. The inclusion of lensing has a dramatic 
effect (by orders of magnitude) on the expected abundance of such hosts at early cosmic 
times. The prominence and duty cycle of such hosts has important implications^'^ with 
respect to the question of whether the cosmic neutral hydrogen, a cold remnant from the 
big bang, was re- ionized by star light or by quasars^ 

To calculate the probability of gravitational lensing in a quasar sample we need to 
specify the luminosity function (number per comoving volume per unit luminosity) that 
describes the quasar population. At 2; ^ 3 this function is well described by a double 
power-law whose shape does not evolve with redshift^^. 



The observed evolution of the break luminosity is described by the dependence^ 



We find that an intrinsic luminosity function having the parameters (f)^ = 624 Gpc 



Pi = 1.64, Ph = 3.43, L^(0) = 1.5 x IO^^Lq, ag = -0.5, = 1.45, ^ = 2.9, C = 2.7, 
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and the inclusion of gravitational lensing (described below) adequately describes three 
observables, namely the luminosity function at 2; ^ 3, and the number density of quasars 
with absolute B-magnitude Mb < —26 at 2; ~ 4.3 (measured by SDSS from a catalog of 39 
quasars^^ with a median redshift of 2; ~ 4.3) and Mb < —27.6 at 2; ~ 6.0. The parameter 
aq is the slope assumed for the typical quasar continuum L(i/) oc L"'". We are interested in 
the number of quasars with luminosities higher than the limiting magnitude z^i^^^ = 20.2, 
which is A^(> I/um, z) — dL (j){L, z) where Li^a is the luminosity of a quasar at redshift 
z corresponding to an apparent magnitude z*:^^. was determined from z*^^ using a 
luminosity distance and a fc-correction computed from a model quasar spectrum including 
the mean absorption by the intergalactic medium^^. 

Gravitational lensing is expected to be highly probable for very luminous quasars^^. 
To illustrate this point, we consider a fictitious gravitational lens that always produces a 
magnification of // = 4 for the sum of multiple images [the average value for a singular 
isothermal sphere (SIS)] but ji — 1 otherwise. We define Tmuit as the probability that a 
random quasar selected in the source plane will be multiply imaged^^, and Fmi to be 
the magnification biased probability that an observed quasar will be multiply imaged. 
Surveys for quasars at 2; < 3 have limiting magnitudes fainter than the break magnitude 
{niB ~ 19). While Tmuu ~ 0.002 at 2; ~ 2, a survey for quasars to a limit will find a 
number of lensed sources that is larger by a bias factor of A^(< Lu^/ /i, z)/N{< L^^, z). At 
2; = 2, Liiin(2;) <^ L^{z) and for /3; = 1.6 the magnification bias factor is ~ 2.3, resulting in 
Fmi ~ 0.005. At 2; = 6, Tmuu ~ 0.008 is significantly higher^^. Furthermore, the limiting 
magnitude of the z ^ 5.8 survey is significantly brighter than the break magnitude and so 
— 3.4 and the bias factor rises to ~ 28. Under these circumstances, Fmi rises to ~ 0.22. 
These simple arguments are consistent with previous estimates of the lensing rate at high 
redshift and demonstrate that lensing has a strong effect on observations of the bright 
SDSS quasars at 2; ;^ 5.8. 
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To find the magnification bias more accurately, we have computed the probabihty 
distributions ^^^p and for the magnification fj, of randomly positioned singly-imaged 

sources, and for the sum of magnifications /^tot of randomly positioned multiply-imaged 
sources due to gravitational lensing by foreground galaxies. We assume that the lens 
galaxies have a constant co-moving density [as the lensing rate for an evolving (Press- 
Schechtcr) population of lenses differs only by <^ 10% from this case^''] and are primarily 
early- type (E/SO) SIS galaxies^^ whose population is described by a Schechter function 
with parameters'^ rii, — 0.27 x 10~^Mpc~^ and ag — —0.5. We assume the Faber- Jackson 
relation (Lg/Lg^) = {ag/a^)'^ where ag is the velocity dispersion of the lens galaxy, with 
cr^ = 220 km see"' and a dark matter velocity dispersion that equals the stellar velocity 
dispersion'^. We ignore dust extinction by the lens galaxy, which mainly arises in the 
much rarer spiral galaxy lenses. Potential lens galaxies must not be detectable in the 
survey data used to select the objects. Galaxies having i* < 22.2 (around 30% of the 
potential lens population) are not considered part of the lens population. To compute i* 
for a galaxy having velocity dispersion a at redshift z, we use Lg^ from the 2dF early-type 
galaxy luminosity function'^, the Faber- Jackson relation, color transformations with a 
/c-correction^'^'^', and the evolution of the rest- frame B-band mass-to- light ratio^^. 

Our lens model includes microlensing by the population of galactic stars which is 
modeled as a de Vaucouleurs profile of point masses embedded in the overall SIS mass 
distribution. The surface mass density in stars for galaxies at 2; = is normalized so 
that the total cosmological density parameter of stars^^ equals 0.005. At 2; > the total 
mass density in stars is assumed to be proportional to the cumulative star-formation 
history^^'^^. The parameters of the de Vaucouleurs profiles are taken from a study of the 
fundamental plane^^, the microlens mass is chosen as O.IM© and the source size as lO'^cm 
(corresponding to ten Schwarzschild radii of a 3 x IO^Mq black hole). The probabilities 
"^-^i^ and ^imii closely resemble the standard form for the SIS and were computed^^ by 
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combining numerical magnification maps with the distribution of microlensing optical depth 



and shear along lensed fines of sight, and the distribution 
was normalized to have unit mean. 



Tmuh J,, ' \^ ^mult 



dtx 'V- 'vault J dfj. 



The fraction of sources which are multiply imaged due to gravitational lensing is 

c dn' + (1 - w)^] A^(> ^,zy 

where Tmuit = 0.0059. We find a value of Fmi ~ 0.30 for the color-selected, fiux-fimited 
z ^ 5.8 sample. This value is higher by two orders of magnitude than the lens fraction at 
low redshifts and demonstrates that lensing must already be considered in a sample with 
only 4 objects. For comparison, our model (including the entire lens population of lens 
galaxies) predicts Fmi ~ 0.01 at ;i; ~ 2 for ttib < 20. These calculations do not include 
selection effects for flux ratios (though the large magnification bias will favor flux ratios 
near 1) or image separations which may serve to lower Fmi- 

A sub arc-second resolution K-band image has been obtained for the z — 5.80 quasar^, 
and it was found to be an unresolved point source. To our knowledge this is the only 2; ;^ 5.8 
quasar for which sub arc-second resolution optical imaging is currently available. However 
a program to image these quasars with the Hubble Space Telescope which will determine 
the multiple-image fraction, is expected to begin within the coming year (X. Fan private 
communication) . Note though that single image quasars might still be magnified by a factor 
of ~ 2. Recent Chandra observations^^ of the z = 6.28 quasar show photons detected on 
the edge of the extraction (1.2") circle, offering a hint that this quasar may be lensed. 

We have computed the distribution of magnifications observed for a sample of quasars 
brighter than Liim at redshift z, 
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In Figure 2 we show the probabihty that the magnification of a quasar is higher than //obs) 
assuming that the quasar belongs to a sample at a redshift 2; ~ 6 with the SDSS magnitude 
limit of z* < 20.2. The plotted distribution is highly skewed; the median magnification is 
med{^obs) ~ 1-2 while the mean magnification is as high as (/x) = 24. Thus, one or more of 
the 2; ;^ 5.8 quasars are likely to be highly magnified while most should be magnified at a 
low level. 

The large values of the magnifications and the highly skewed shape of the distributions 
in Figure 2 suggest that lensing must alter the observed luminosity function. Indeed, we 
find that the space density of quasars with Mb < —27.6 is increased by 40%, and that the 
slope is decreased by 0.15. A quasar magnified by /i is detectable to luminosities as low 
as Liim/ jJi and has its luminosity, L, overestimated by jj,. The factor Rl^ hj which the 
luminosity density of quasars brighter than is overestimated due to lensing is therefore 



We find Rld ~ 2, implying that naive computation of the quasar luminosity density from 
the 2; ;^ 5.8 sample might significantly overestimate its true value. Magnification bias also 
affects quasars that are not multiply imaged. We therefore predict an enhanced angular 
correlation on the sky between 2; ~ 6 quasars and foreground galaxies. 

The high source redshifts imply image separations that are shghtly larger than usuaP^, 
~ 1-2". Furthermore, the lenses are likely to be found at higher redshift^^ due in part to 
the higher source redshifts but also because bright (low- 2;) lenses are excluded. Lensing of 
high-redshift quasars allows measurement of the masses of the lens galaxies^^ at redshifts 
higher than currently possible. Furthermore, quasar microlensing should be common over 
a ten year baseline^''. This offers the exciting possibility of measuring the source size, 
and hence black-hole mass, which in turn yields the ratio between quasar luminosity and 
its Eddington value, a quantity that will be useful in constraining models of structure 
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formation^. 

The recently published spectra of the highest redshift quasar ^'^^ limit the flux in the 
Gunn-Peterson^° trough to less than 3 x 10~^^ ergsec"^ A~^. Given the importance of this 
first observation of the re-ionization epoch and of future complimentary observations, we are 
motivated to ask whether or not lens galaxy light will contaminate deep observations of the 
Gunn-Peterson trough in the quasar spectrum, even though such a galaxy is not detected 
in the initial imaging survey. We have computed the flux of lens galaxies at redshift z with 
velocity dispersion o", and convolved the results with the joint probability distribution for 
the lens galaxy redshift and velocity dispersion. We find that ~ 40% of multiple image 
lens galaxies {%* < 22.2) will contribute flux in the Gunn-Peterson trough above a level of 
3 X 10~^^ ergsec"^ A~^. For some quasars the contamination of the Gunn-Peterson trough 
by flux from lens galaxies may limit the ability of deep spectroscopic observations to probe 
the evolution of the neutral hydrogen fraction during the epoch of reionization. 

Correspondence and requests for materials to Abraham Loeb. 
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Fig. 1. — Enhancement in the co- moving space density of host galaxies for quasars at 2; = 6 
as a function of the magnification due to lensing, fj,. Ignoring lensing, the Eddington hmit 
imphes^ a minimum black hole mass of ~ 3 x lO^M© for the SDSS quasars at 2 ~ 6. 
The inclusion of lensing reduces the implied (minimum) black hole mass by a factor of /i. 
Given some fixed efficiency for assembling gas into a central black hole within a galaxy, the 
implied mass of the host galaxy is also lowered by n after lensing is included. Plotted is 
^(/^) = d{iogW/^^\) ^di^v wliere ^^1^ is the Press- Schechter^ comoving density of galaxy 
halos with mass M per logarithmic interval in M. We show curves for three possible halo 
masses of quasar hosts, Mhaio = lO^^M©, 10^'^ M©, and lO^^M©. As seen, the enhancement in 
the implied abundance of quasar hosts increases with increasing halo mass. The above three 
choices for halo masses are all conservatively lower than inferred for the same black hole 
mass in the local universe^. Throughout this Letter we assume the standard cosmological 
parameters = 0.35, = 0.65, ^2;, = 0.05, Hq = 65kmsec~^ Mpc, n = 1 and as = 0.8. 
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Fig. 2. — The probability of observing a magnification larger than //obg for a quasar at a 
redshift z = 6 in a sample with a magnitude limit z* < 20.2. The distribution is highly 
skewed, having a median of med{nobs) = 1-2 and a mean of (fiobs) = 24.0. The multiple 
image fraction is Fmi — 0.30. We have also computed the a-posteriori values of Fmi and 
(//) for specific quasars. For SDSS 0836-0054 {z = 5.82), SDSS 1306-0356 {z = 5.99) and 
SDSS 1030-0524 {z = 6.28) we find Fmi = 0.40, 0.32 and 0.31, and (fi) = 50, 25 and 23, 
respectively. While Fan et al.^ find that (3^ = 3.43 is consistent with the luminosities of 
the z ^ 5.8 quasars, it is possible that the luminosity function at high redshift is not as 
steep as /3h — 3.43. In this case, the magnification bias will not be as large and the mean 
magnification and fraction of quasars that are multiply imaged will be lower. Wc have re- 
computed the lens statistics assuming that the bright end slope is significantly fiatter at high 
redshift. Assuming that Pi — 1.64 at all redshifts and Ph — 3.43 for 2; < 3 but Ph — 2.58 for 
z > 3 (the value found^^ for quasars at 2; ~ 4.3) we inferred similar parameters to describe 
the observed luminosity function as before [I'*(0) = 1.5 x lO^^L©, z^, = 1.6, ^ = 3.3, ( = 2.65]. 
Remarkably, the multiple image fraction is still nearly 0.1 in this case. Additional uncertainty 
in the calculation arises in the choice of lens model. The value Tmuit is proportional to n^a^. 
The dependence of Fmi on Tmuit is complex [see equation (3)], however large magnification 
biases result in a relation that is less sensitive than linear. For example, reducing the value 
of Tmuit by a factor of 1.5 results in Fmi = 0.22 if Ph = 3.43 and Fmi = 0.043 if Ph = 2.58. 



